\(\int \text {csch}^3(c+d x) (a+b \tanh ^3(c+d x))^2 \, dx\) [63]
Optimal result
Integrand size = 23, antiderivative size = 107 \[
\int \text {csch}^3(c+d x) \left (a+b \tanh ^3(c+d x)\right )^2 \, dx=\frac {a b \arctan (\sinh (c+d x))}{d}+\frac {a^2 \text {arctanh}(\cosh (c+d x))}{2 d}-\frac {a^2 \coth (c+d x) \text {csch}(c+d x)}{2 d}-\frac {b^2 \text {sech}^3(c+d x)}{3 d}+\frac {b^2 \text {sech}^5(c+d x)}{5 d}+\frac {a b \text {sech}(c+d x) \tanh (c+d x)}{d}
\]
[Out]
a*b*arctan(sinh(d*x+c))/d+1/2*a^2*arctanh(cosh(d*x+c))/d-1/2*a^2*coth(d*x+c)*csch(d*x+c)/d-1/3*b^2*sech(d*x+c)
^3/d+1/5*b^2*sech(d*x+c)^5/d+a*b*sech(d*x+c)*tanh(d*x+c)/d
Rubi [A] (verified)
Time = 0.12 (sec) , antiderivative size = 107, normalized size of antiderivative = 1.00, number of
steps used = 9, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.217, Rules used = {3747, 3853, 3855, 2686, 14}
\[
\int \text {csch}^3(c+d x) \left (a+b \tanh ^3(c+d x)\right )^2 \, dx=\frac {a^2 \text {arctanh}(\cosh (c+d x))}{2 d}-\frac {a^2 \coth (c+d x) \text {csch}(c+d x)}{2 d}+\frac {a b \arctan (\sinh (c+d x))}{d}+\frac {a b \tanh (c+d x) \text {sech}(c+d x)}{d}+\frac {b^2 \text {sech}^5(c+d x)}{5 d}-\frac {b^2 \text {sech}^3(c+d x)}{3 d}
\]
[In]
Int[Csch[c + d*x]^3*(a + b*Tanh[c + d*x]^3)^2,x]
[Out]
(a*b*ArcTan[Sinh[c + d*x]])/d + (a^2*ArcTanh[Cosh[c + d*x]])/(2*d) - (a^2*Coth[c + d*x]*Csch[c + d*x])/(2*d) -
(b^2*Sech[c + d*x]^3)/(3*d) + (b^2*Sech[c + d*x]^5)/(5*d) + (a*b*Sech[c + d*x]*Tanh[c + d*x])/d
Rule 14
Int[(u_)*((c_.)*(x_))^(m_.), x_Symbol] :> Int[ExpandIntegrand[(c*x)^m*u, x], x] /; FreeQ[{c, m}, x] && SumQ[u]
&& !LinearQ[u, x] && !MatchQ[u, (a_) + (b_.)*(v_) /; FreeQ[{a, b}, x] && InverseFunctionQ[v]]
Rule 2686
Int[((a_.)*sec[(e_.) + (f_.)*(x_)])^(m_.)*((b_.)*tan[(e_.) + (f_.)*(x_)])^(n_.), x_Symbol] :> Dist[a/f, Subst[
Int[(a*x)^(m - 1)*(-1 + x^2)^((n - 1)/2), x], x, Sec[e + f*x]], x] /; FreeQ[{a, e, f, m}, x] && IntegerQ[(n -
1)/2] && !(IntegerQ[m/2] && LtQ[0, m, n + 1])
Rule 3747
Int[((d_.)*sin[(e_.) + (f_.)*(x_)])^(m_.)*((a_) + (b_.)*((c_.)*tan[(e_.) + (f_.)*(x_)])^(n_))^(p_.), x_Symbol]
:> Int[ExpandTrig[(d*sin[e + f*x])^m*(a + b*(c*tan[e + f*x])^n)^p, x], x] /; FreeQ[{a, b, c, d, e, f, m, n},
x] && IGtQ[p, 0]
Rule 3853
Int[(csc[(c_.) + (d_.)*(x_)]*(b_.))^(n_), x_Symbol] :> Simp[(-b)*Cos[c + d*x]*((b*Csc[c + d*x])^(n - 1)/(d*(n
- 1))), x] + Dist[b^2*((n - 2)/(n - 1)), Int[(b*Csc[c + d*x])^(n - 2), x], x] /; FreeQ[{b, c, d}, x] && GtQ[n,
1] && IntegerQ[2*n]
Rule 3855
Int[csc[(c_.) + (d_.)*(x_)], x_Symbol] :> Simp[-ArcTanh[Cos[c + d*x]]/d, x] /; FreeQ[{c, d}, x]
Rubi steps \begin{align*}
\text {integral}& = -\left (i \int \left (i a^2 \text {csch}^3(c+d x)+2 i a b \text {sech}^3(c+d x)+i b^2 \text {sech}^3(c+d x) \tanh ^3(c+d x)\right ) \, dx\right ) \\ & = a^2 \int \text {csch}^3(c+d x) \, dx+(2 a b) \int \text {sech}^3(c+d x) \, dx+b^2 \int \text {sech}^3(c+d x) \tanh ^3(c+d x) \, dx \\ & = -\frac {a^2 \coth (c+d x) \text {csch}(c+d x)}{2 d}+\frac {a b \text {sech}(c+d x) \tanh (c+d x)}{d}-\frac {1}{2} a^2 \int \text {csch}(c+d x) \, dx+(a b) \int \text {sech}(c+d x) \, dx+\frac {b^2 \text {Subst}\left (\int x^2 \left (-1+x^2\right ) \, dx,x,\text {sech}(c+d x)\right )}{d} \\ & = \frac {a b \arctan (\sinh (c+d x))}{d}+\frac {a^2 \text {arctanh}(\cosh (c+d x))}{2 d}-\frac {a^2 \coth (c+d x) \text {csch}(c+d x)}{2 d}+\frac {a b \text {sech}(c+d x) \tanh (c+d x)}{d}+\frac {b^2 \text {Subst}\left (\int \left (-x^2+x^4\right ) \, dx,x,\text {sech}(c+d x)\right )}{d} \\ & = \frac {a b \arctan (\sinh (c+d x))}{d}+\frac {a^2 \text {arctanh}(\cosh (c+d x))}{2 d}-\frac {a^2 \coth (c+d x) \text {csch}(c+d x)}{2 d}-\frac {b^2 \text {sech}^3(c+d x)}{3 d}+\frac {b^2 \text {sech}^5(c+d x)}{5 d}+\frac {a b \text {sech}(c+d x) \tanh (c+d x)}{d} \\
\end{align*}
Mathematica [A] (verified)
Time = 1.74 (sec) , antiderivative size = 159, normalized size of antiderivative = 1.49
\[
\int \text {csch}^3(c+d x) \left (a+b \tanh ^3(c+d x)\right )^2 \, dx=\frac {2 a b \arctan \left (\tanh \left (\frac {1}{2} (c+d x)\right )\right )}{d}-\frac {a^2 \text {csch}^2\left (\frac {1}{2} (c+d x)\right )}{8 d}+\frac {a^2 \log \left (\cosh \left (\frac {1}{2} (c+d x)\right )\right )}{2 d}-\frac {a^2 \log \left (\sinh \left (\frac {1}{2} (c+d x)\right )\right )}{2 d}-\frac {a^2 \text {sech}^2\left (\frac {1}{2} (c+d x)\right )}{8 d}-\frac {b^2 \text {sech}^3(c+d x)}{3 d}+\frac {b^2 \text {sech}^5(c+d x)}{5 d}+\frac {a b \text {sech}(c+d x) \tanh (c+d x)}{d}
\]
[In]
Integrate[Csch[c + d*x]^3*(a + b*Tanh[c + d*x]^3)^2,x]
[Out]
(2*a*b*ArcTan[Tanh[(c + d*x)/2]])/d - (a^2*Csch[(c + d*x)/2]^2)/(8*d) + (a^2*Log[Cosh[(c + d*x)/2]])/(2*d) - (
a^2*Log[Sinh[(c + d*x)/2]])/(2*d) - (a^2*Sech[(c + d*x)/2]^2)/(8*d) - (b^2*Sech[c + d*x]^3)/(3*d) + (b^2*Sech[
c + d*x]^5)/(5*d) + (a*b*Sech[c + d*x]*Tanh[c + d*x])/d
Maple [A] (verified)
Time = 8.68 (sec) , antiderivative size = 91, normalized size of antiderivative = 0.85
| | |
| method | result | size |
| | |
| derivativedivides |
\(\frac {a^{2} \left (-\frac {\operatorname {csch}\left (d x +c \right ) \coth \left (d x +c \right )}{2}+\operatorname {arctanh}\left ({\mathrm e}^{d x +c}\right )\right )+2 a b \left (\frac {\operatorname {sech}\left (d x +c \right ) \tanh \left (d x +c \right )}{2}+\arctan \left ({\mathrm e}^{d x +c}\right )\right )+b^{2} \left (-\frac {\sinh \left (d x +c \right )^{2}}{3 \cosh \left (d x +c \right )^{5}}-\frac {2}{15 \cosh \left (d x +c \right )^{5}}\right )}{d}\) |
\(91\) |
| default |
\(\frac {a^{2} \left (-\frac {\operatorname {csch}\left (d x +c \right ) \coth \left (d x +c \right )}{2}+\operatorname {arctanh}\left ({\mathrm e}^{d x +c}\right )\right )+2 a b \left (\frac {\operatorname {sech}\left (d x +c \right ) \tanh \left (d x +c \right )}{2}+\arctan \left ({\mathrm e}^{d x +c}\right )\right )+b^{2} \left (-\frac {\sinh \left (d x +c \right )^{2}}{3 \cosh \left (d x +c \right )^{5}}-\frac {2}{15 \cosh \left (d x +c \right )^{5}}\right )}{d}\) |
\(91\) |
| risch |
\(-\frac {{\mathrm e}^{d x +c} \left (15 a^{2} {\mathrm e}^{12 d x +12 c}-30 a b \,{\mathrm e}^{12 d x +12 c}+90 a^{2} {\mathrm e}^{10 d x +10 c}+40 b^{2} {\mathrm e}^{10 d x +10 c}+225 a^{2} {\mathrm e}^{8 d x +8 c}+90 a b \,{\mathrm e}^{8 d x +8 c}-96 b^{2} {\mathrm e}^{8 d x +8 c}+300 a^{2} {\mathrm e}^{6 d x +6 c}+112 b^{2} {\mathrm e}^{6 d x +6 c}+225 a^{2} {\mathrm e}^{4 d x +4 c}-90 a b \,{\mathrm e}^{4 d x +4 c}-96 \,{\mathrm e}^{4 d x +4 c} b^{2}+90 a^{2} {\mathrm e}^{2 d x +2 c}+40 \,{\mathrm e}^{2 d x +2 c} b^{2}+15 a^{2}+30 a b \right )}{15 d \left ({\mathrm e}^{2 d x +2 c}-1\right )^{2} \left ({\mathrm e}^{2 d x +2 c}+1\right )^{5}}+\frac {a^{2} \ln \left ({\mathrm e}^{d x +c}+1\right )}{2 d}-\frac {a^{2} \ln \left ({\mathrm e}^{d x +c}-1\right )}{2 d}+\frac {i b a \ln \left ({\mathrm e}^{d x +c}+i\right )}{d}-\frac {i b a \ln \left ({\mathrm e}^{d x +c}-i\right )}{d}\) |
\(312\) |
| | |
|
|
|
|
[In]
int(csch(d*x+c)^3*(a+b*tanh(d*x+c)^3)^2,x,method=_RETURNVERBOSE)
[Out]
1/d*(a^2*(-1/2*csch(d*x+c)*coth(d*x+c)+arctanh(exp(d*x+c)))+2*a*b*(1/2*sech(d*x+c)*tanh(d*x+c)+arctan(exp(d*x+
c)))+b^2*(-1/3*sinh(d*x+c)^2/cosh(d*x+c)^5-2/15/cosh(d*x+c)^5))
Fricas [B] (verification not implemented)
Leaf count of result is larger than twice the leaf count of optimal. 4642 vs. \(2 (99) = 198\).
Time = 0.32 (sec) , antiderivative size = 4642, normalized size of antiderivative = 43.38
\[
\int \text {csch}^3(c+d x) \left (a+b \tanh ^3(c+d x)\right )^2 \, dx=\text {Too large to display}
\]
[In]
integrate(csch(d*x+c)^3*(a+b*tanh(d*x+c)^3)^2,x, algorithm="fricas")
[Out]
-1/30*(30*(a^2 - 2*a*b)*cosh(d*x + c)^13 + 390*(a^2 - 2*a*b)*cosh(d*x + c)*sinh(d*x + c)^12 + 30*(a^2 - 2*a*b)
*sinh(d*x + c)^13 + 20*(9*a^2 + 4*b^2)*cosh(d*x + c)^11 + 20*(117*(a^2 - 2*a*b)*cosh(d*x + c)^2 + 9*a^2 + 4*b^
2)*sinh(d*x + c)^11 + 220*(39*(a^2 - 2*a*b)*cosh(d*x + c)^3 + (9*a^2 + 4*b^2)*cosh(d*x + c))*sinh(d*x + c)^10
+ 6*(75*a^2 + 30*a*b - 32*b^2)*cosh(d*x + c)^9 + 2*(10725*(a^2 - 2*a*b)*cosh(d*x + c)^4 + 550*(9*a^2 + 4*b^2)*
cosh(d*x + c)^2 + 225*a^2 + 90*a*b - 96*b^2)*sinh(d*x + c)^9 + 6*(6435*(a^2 - 2*a*b)*cosh(d*x + c)^5 + 550*(9*
a^2 + 4*b^2)*cosh(d*x + c)^3 + 9*(75*a^2 + 30*a*b - 32*b^2)*cosh(d*x + c))*sinh(d*x + c)^8 + 8*(75*a^2 + 28*b^
2)*cosh(d*x + c)^7 + 8*(6435*(a^2 - 2*a*b)*cosh(d*x + c)^6 + 825*(9*a^2 + 4*b^2)*cosh(d*x + c)^4 + 27*(75*a^2
+ 30*a*b - 32*b^2)*cosh(d*x + c)^2 + 75*a^2 + 28*b^2)*sinh(d*x + c)^7 + 8*(6435*(a^2 - 2*a*b)*cosh(d*x + c)^7
+ 1155*(9*a^2 + 4*b^2)*cosh(d*x + c)^5 + 63*(75*a^2 + 30*a*b - 32*b^2)*cosh(d*x + c)^3 + 7*(75*a^2 + 28*b^2)*c
osh(d*x + c))*sinh(d*x + c)^6 + 6*(75*a^2 - 30*a*b - 32*b^2)*cosh(d*x + c)^5 + 6*(6435*(a^2 - 2*a*b)*cosh(d*x
+ c)^8 + 1540*(9*a^2 + 4*b^2)*cosh(d*x + c)^6 + 126*(75*a^2 + 30*a*b - 32*b^2)*cosh(d*x + c)^4 + 28*(75*a^2 +
28*b^2)*cosh(d*x + c)^2 + 75*a^2 - 30*a*b - 32*b^2)*sinh(d*x + c)^5 + 2*(10725*(a^2 - 2*a*b)*cosh(d*x + c)^9 +
3300*(9*a^2 + 4*b^2)*cosh(d*x + c)^7 + 378*(75*a^2 + 30*a*b - 32*b^2)*cosh(d*x + c)^5 + 140*(75*a^2 + 28*b^2)
*cosh(d*x + c)^3 + 15*(75*a^2 - 30*a*b - 32*b^2)*cosh(d*x + c))*sinh(d*x + c)^4 + 20*(9*a^2 + 4*b^2)*cosh(d*x
+ c)^3 + 4*(2145*(a^2 - 2*a*b)*cosh(d*x + c)^10 + 825*(9*a^2 + 4*b^2)*cosh(d*x + c)^8 + 126*(75*a^2 + 30*a*b -
32*b^2)*cosh(d*x + c)^6 + 70*(75*a^2 + 28*b^2)*cosh(d*x + c)^4 + 15*(75*a^2 - 30*a*b - 32*b^2)*cosh(d*x + c)^
2 + 45*a^2 + 20*b^2)*sinh(d*x + c)^3 + 4*(585*(a^2 - 2*a*b)*cosh(d*x + c)^11 + 275*(9*a^2 + 4*b^2)*cosh(d*x +
c)^9 + 54*(75*a^2 + 30*a*b - 32*b^2)*cosh(d*x + c)^7 + 42*(75*a^2 + 28*b^2)*cosh(d*x + c)^5 + 15*(75*a^2 - 30*
a*b - 32*b^2)*cosh(d*x + c)^3 + 15*(9*a^2 + 4*b^2)*cosh(d*x + c))*sinh(d*x + c)^2 - 60*(a*b*cosh(d*x + c)^14 +
14*a*b*cosh(d*x + c)*sinh(d*x + c)^13 + a*b*sinh(d*x + c)^14 + 3*a*b*cosh(d*x + c)^12 + (91*a*b*cosh(d*x + c)
^2 + 3*a*b)*sinh(d*x + c)^12 + a*b*cosh(d*x + c)^10 + 4*(91*a*b*cosh(d*x + c)^3 + 9*a*b*cosh(d*x + c))*sinh(d*
x + c)^11 + (1001*a*b*cosh(d*x + c)^4 + 198*a*b*cosh(d*x + c)^2 + a*b)*sinh(d*x + c)^10 - 5*a*b*cosh(d*x + c)^
8 + 2*(1001*a*b*cosh(d*x + c)^5 + 330*a*b*cosh(d*x + c)^3 + 5*a*b*cosh(d*x + c))*sinh(d*x + c)^9 + (3003*a*b*c
osh(d*x + c)^6 + 1485*a*b*cosh(d*x + c)^4 + 45*a*b*cosh(d*x + c)^2 - 5*a*b)*sinh(d*x + c)^8 - 5*a*b*cosh(d*x +
c)^6 + 8*(429*a*b*cosh(d*x + c)^7 + 297*a*b*cosh(d*x + c)^5 + 15*a*b*cosh(d*x + c)^3 - 5*a*b*cosh(d*x + c))*s
inh(d*x + c)^7 + (3003*a*b*cosh(d*x + c)^8 + 2772*a*b*cosh(d*x + c)^6 + 210*a*b*cosh(d*x + c)^4 - 140*a*b*cosh
(d*x + c)^2 - 5*a*b)*sinh(d*x + c)^6 + a*b*cosh(d*x + c)^4 + 2*(1001*a*b*cosh(d*x + c)^9 + 1188*a*b*cosh(d*x +
c)^7 + 126*a*b*cosh(d*x + c)^5 - 140*a*b*cosh(d*x + c)^3 - 15*a*b*cosh(d*x + c))*sinh(d*x + c)^5 + (1001*a*b*
cosh(d*x + c)^10 + 1485*a*b*cosh(d*x + c)^8 + 210*a*b*cosh(d*x + c)^6 - 350*a*b*cosh(d*x + c)^4 - 75*a*b*cosh(
d*x + c)^2 + a*b)*sinh(d*x + c)^4 + 3*a*b*cosh(d*x + c)^2 + 4*(91*a*b*cosh(d*x + c)^11 + 165*a*b*cosh(d*x + c)
^9 + 30*a*b*cosh(d*x + c)^7 - 70*a*b*cosh(d*x + c)^5 - 25*a*b*cosh(d*x + c)^3 + a*b*cosh(d*x + c))*sinh(d*x +
c)^3 + (91*a*b*cosh(d*x + c)^12 + 198*a*b*cosh(d*x + c)^10 + 45*a*b*cosh(d*x + c)^8 - 140*a*b*cosh(d*x + c)^6
- 75*a*b*cosh(d*x + c)^4 + 6*a*b*cosh(d*x + c)^2 + 3*a*b)*sinh(d*x + c)^2 + a*b + 2*(7*a*b*cosh(d*x + c)^13 +
18*a*b*cosh(d*x + c)^11 + 5*a*b*cosh(d*x + c)^9 - 20*a*b*cosh(d*x + c)^7 - 15*a*b*cosh(d*x + c)^5 + 2*a*b*cosh
(d*x + c)^3 + 3*a*b*cosh(d*x + c))*sinh(d*x + c))*arctan(cosh(d*x + c) + sinh(d*x + c)) + 30*(a^2 + 2*a*b)*cos
h(d*x + c) - 15*(a^2*cosh(d*x + c)^14 + 14*a^2*cosh(d*x + c)*sinh(d*x + c)^13 + a^2*sinh(d*x + c)^14 + 3*a^2*c
osh(d*x + c)^12 + (91*a^2*cosh(d*x + c)^2 + 3*a^2)*sinh(d*x + c)^12 + a^2*cosh(d*x + c)^10 + 4*(91*a^2*cosh(d*
x + c)^3 + 9*a^2*cosh(d*x + c))*sinh(d*x + c)^11 + (1001*a^2*cosh(d*x + c)^4 + 198*a^2*cosh(d*x + c)^2 + a^2)*
sinh(d*x + c)^10 - 5*a^2*cosh(d*x + c)^8 + 2*(1001*a^2*cosh(d*x + c)^5 + 330*a^2*cosh(d*x + c)^3 + 5*a^2*cosh(
d*x + c))*sinh(d*x + c)^9 + (3003*a^2*cosh(d*x + c)^6 + 1485*a^2*cosh(d*x + c)^4 + 45*a^2*cosh(d*x + c)^2 - 5*
a^2)*sinh(d*x + c)^8 - 5*a^2*cosh(d*x + c)^6 + 8*(429*a^2*cosh(d*x + c)^7 + 297*a^2*cosh(d*x + c)^5 + 15*a^2*c
osh(d*x + c)^3 - 5*a^2*cosh(d*x + c))*sinh(d*x + c)^7 + (3003*a^2*cosh(d*x + c)^8 + 2772*a^2*cosh(d*x + c)^6 +
210*a^2*cosh(d*x + c)^4 - 140*a^2*cosh(d*x + c)^2 - 5*a^2)*sinh(d*x + c)^6 + a^2*cosh(d*x + c)^4 + 2*(1001*a^
2*cosh(d*x + c)^9 + 1188*a^2*cosh(d*x + c)^7 + 126*a^2*cosh(d*x + c)^5 - 140*a^2*cosh(d*x + c)^3 - 15*a^2*cosh
(d*x + c))*sinh(d*x + c)^5 + (1001*a^2*cosh(d*x + c)^10 + 1485*a^2*cosh(d*x + c)^8 + 210*a^2*cosh(d*x + c)^6 -
350*a^2*cosh(d*x + c)^4 - 75*a^2*cosh(d*x + c)^2 + a^2)*sinh(d*x + c)^4 + 3*a^2*cosh(d*x + c)^2 + 4*(91*a^2*c
osh(d*x + c)^11 + 165*a^2*cosh(d*x + c)^9 + 30*a^2*cosh(d*x + c)^7 - 70*a^2*cosh(d*x + c)^5 - 25*a^2*cosh(d*x
+ c)^3 + a^2*cosh(d*x + c))*sinh(d*x + c)^3 + (91*a^2*cosh(d*x + c)^12 + 198*a^2*cosh(d*x + c)^10 + 45*a^2*cos
h(d*x + c)^8 - 140*a^2*cosh(d*x + c)^6 - 75*a^2*cosh(d*x + c)^4 + 6*a^2*cosh(d*x + c)^2 + 3*a^2)*sinh(d*x + c)
^2 + a^2 + 2*(7*a^2*cosh(d*x + c)^13 + 18*a^2*cosh(d*x + c)^11 + 5*a^2*cosh(d*x + c)^9 - 20*a^2*cosh(d*x + c)^
7 - 15*a^2*cosh(d*x + c)^5 + 2*a^2*cosh(d*x + c)^3 + 3*a^2*cosh(d*x + c))*sinh(d*x + c))*log(cosh(d*x + c) + s
inh(d*x + c) + 1) + 15*(a^2*cosh(d*x + c)^14 + 14*a^2*cosh(d*x + c)*sinh(d*x + c)^13 + a^2*sinh(d*x + c)^14 +
3*a^2*cosh(d*x + c)^12 + (91*a^2*cosh(d*x + c)^2 + 3*a^2)*sinh(d*x + c)^12 + a^2*cosh(d*x + c)^10 + 4*(91*a^2*
cosh(d*x + c)^3 + 9*a^2*cosh(d*x + c))*sinh(d*x + c)^11 + (1001*a^2*cosh(d*x + c)^4 + 198*a^2*cosh(d*x + c)^2
+ a^2)*sinh(d*x + c)^10 - 5*a^2*cosh(d*x + c)^8 + 2*(1001*a^2*cosh(d*x + c)^5 + 330*a^2*cosh(d*x + c)^3 + 5*a^
2*cosh(d*x + c))*sinh(d*x + c)^9 + (3003*a^2*cosh(d*x + c)^6 + 1485*a^2*cosh(d*x + c)^4 + 45*a^2*cosh(d*x + c)
^2 - 5*a^2)*sinh(d*x + c)^8 - 5*a^2*cosh(d*x + c)^6 + 8*(429*a^2*cosh(d*x + c)^7 + 297*a^2*cosh(d*x + c)^5 + 1
5*a^2*cosh(d*x + c)^3 - 5*a^2*cosh(d*x + c))*sinh(d*x + c)^7 + (3003*a^2*cosh(d*x + c)^8 + 2772*a^2*cosh(d*x +
c)^6 + 210*a^2*cosh(d*x + c)^4 - 140*a^2*cosh(d*x + c)^2 - 5*a^2)*sinh(d*x + c)^6 + a^2*cosh(d*x + c)^4 + 2*(
1001*a^2*cosh(d*x + c)^9 + 1188*a^2*cosh(d*x + c)^7 + 126*a^2*cosh(d*x + c)^5 - 140*a^2*cosh(d*x + c)^3 - 15*a
^2*cosh(d*x + c))*sinh(d*x + c)^5 + (1001*a^2*cosh(d*x + c)^10 + 1485*a^2*cosh(d*x + c)^8 + 210*a^2*cosh(d*x +
c)^6 - 350*a^2*cosh(d*x + c)^4 - 75*a^2*cosh(d*x + c)^2 + a^2)*sinh(d*x + c)^4 + 3*a^2*cosh(d*x + c)^2 + 4*(9
1*a^2*cosh(d*x + c)^11 + 165*a^2*cosh(d*x + c)^9 + 30*a^2*cosh(d*x + c)^7 - 70*a^2*cosh(d*x + c)^5 - 25*a^2*co
sh(d*x + c)^3 + a^2*cosh(d*x + c))*sinh(d*x + c)^3 + (91*a^2*cosh(d*x + c)^12 + 198*a^2*cosh(d*x + c)^10 + 45*
a^2*cosh(d*x + c)^8 - 140*a^2*cosh(d*x + c)^6 - 75*a^2*cosh(d*x + c)^4 + 6*a^2*cosh(d*x + c)^2 + 3*a^2)*sinh(d
*x + c)^2 + a^2 + 2*(7*a^2*cosh(d*x + c)^13 + 18*a^2*cosh(d*x + c)^11 + 5*a^2*cosh(d*x + c)^9 - 20*a^2*cosh(d*
x + c)^7 - 15*a^2*cosh(d*x + c)^5 + 2*a^2*cosh(d*x + c)^3 + 3*a^2*cosh(d*x + c))*sinh(d*x + c))*log(cosh(d*x +
c) + sinh(d*x + c) - 1) + 2*(195*(a^2 - 2*a*b)*cosh(d*x + c)^12 + 110*(9*a^2 + 4*b^2)*cosh(d*x + c)^10 + 27*(
75*a^2 + 30*a*b - 32*b^2)*cosh(d*x + c)^8 + 28*(75*a^2 + 28*b^2)*cosh(d*x + c)^6 + 15*(75*a^2 - 30*a*b - 32*b^
2)*cosh(d*x + c)^4 + 30*(9*a^2 + 4*b^2)*cosh(d*x + c)^2 + 15*a^2 + 30*a*b)*sinh(d*x + c))/(d*cosh(d*x + c)^14
+ 14*d*cosh(d*x + c)*sinh(d*x + c)^13 + d*sinh(d*x + c)^14 + 3*d*cosh(d*x + c)^12 + (91*d*cosh(d*x + c)^2 + 3*
d)*sinh(d*x + c)^12 + 4*(91*d*cosh(d*x + c)^3 + 9*d*cosh(d*x + c))*sinh(d*x + c)^11 + d*cosh(d*x + c)^10 + (10
01*d*cosh(d*x + c)^4 + 198*d*cosh(d*x + c)^2 + d)*sinh(d*x + c)^10 + 2*(1001*d*cosh(d*x + c)^5 + 330*d*cosh(d*
x + c)^3 + 5*d*cosh(d*x + c))*sinh(d*x + c)^9 - 5*d*cosh(d*x + c)^8 + (3003*d*cosh(d*x + c)^6 + 1485*d*cosh(d*
x + c)^4 + 45*d*cosh(d*x + c)^2 - 5*d)*sinh(d*x + c)^8 + 8*(429*d*cosh(d*x + c)^7 + 297*d*cosh(d*x + c)^5 + 15
*d*cosh(d*x + c)^3 - 5*d*cosh(d*x + c))*sinh(d*x + c)^7 - 5*d*cosh(d*x + c)^6 + (3003*d*cosh(d*x + c)^8 + 2772
*d*cosh(d*x + c)^6 + 210*d*cosh(d*x + c)^4 - 140*d*cosh(d*x + c)^2 - 5*d)*sinh(d*x + c)^6 + 2*(1001*d*cosh(d*x
+ c)^9 + 1188*d*cosh(d*x + c)^7 + 126*d*cosh(d*x + c)^5 - 140*d*cosh(d*x + c)^3 - 15*d*cosh(d*x + c))*sinh(d*
x + c)^5 + d*cosh(d*x + c)^4 + (1001*d*cosh(d*x + c)^10 + 1485*d*cosh(d*x + c)^8 + 210*d*cosh(d*x + c)^6 - 350
*d*cosh(d*x + c)^4 - 75*d*cosh(d*x + c)^2 + d)*sinh(d*x + c)^4 + 4*(91*d*cosh(d*x + c)^11 + 165*d*cosh(d*x + c
)^9 + 30*d*cosh(d*x + c)^7 - 70*d*cosh(d*x + c)^5 - 25*d*cosh(d*x + c)^3 + d*cosh(d*x + c))*sinh(d*x + c)^3 +
3*d*cosh(d*x + c)^2 + (91*d*cosh(d*x + c)^12 + 198*d*cosh(d*x + c)^10 + 45*d*cosh(d*x + c)^8 - 140*d*cosh(d*x
+ c)^6 - 75*d*cosh(d*x + c)^4 + 6*d*cosh(d*x + c)^2 + 3*d)*sinh(d*x + c)^2 + 2*(7*d*cosh(d*x + c)^13 + 18*d*co
sh(d*x + c)^11 + 5*d*cosh(d*x + c)^9 - 20*d*cosh(d*x + c)^7 - 15*d*cosh(d*x + c)^5 + 2*d*cosh(d*x + c)^3 + 3*d
*cosh(d*x + c))*sinh(d*x + c) + d)
Sympy [F]
\[
\int \text {csch}^3(c+d x) \left (a+b \tanh ^3(c+d x)\right )^2 \, dx=\int \left (a + b \tanh ^{3}{\left (c + d x \right )}\right )^{2} \operatorname {csch}^{3}{\left (c + d x \right )}\, dx
\]
[In]
integrate(csch(d*x+c)**3*(a+b*tanh(d*x+c)**3)**2,x)
[Out]
Integral((a + b*tanh(c + d*x)**3)**2*csch(c + d*x)**3, x)
Maxima [B] (verification not implemented)
Leaf count of result is larger than twice the leaf count of optimal. 378 vs. \(2 (99) = 198\).
Time = 0.29 (sec) , antiderivative size = 378, normalized size of antiderivative = 3.53
\[
\int \text {csch}^3(c+d x) \left (a+b \tanh ^3(c+d x)\right )^2 \, dx=-2 \, a b {\left (\frac {\arctan \left (e^{\left (-d x - c\right )}\right )}{d} - \frac {e^{\left (-d x - c\right )} - e^{\left (-3 \, d x - 3 \, c\right )}}{d {\left (2 \, e^{\left (-2 \, d x - 2 \, c\right )} + e^{\left (-4 \, d x - 4 \, c\right )} + 1\right )}}\right )} + \frac {1}{2} \, a^{2} {\left (\frac {\log \left (e^{\left (-d x - c\right )} + 1\right )}{d} - \frac {\log \left (e^{\left (-d x - c\right )} - 1\right )}{d} + \frac {2 \, {\left (e^{\left (-d x - c\right )} + e^{\left (-3 \, d x - 3 \, c\right )}\right )}}{d {\left (2 \, e^{\left (-2 \, d x - 2 \, c\right )} - e^{\left (-4 \, d x - 4 \, c\right )} - 1\right )}}\right )} - \frac {8}{15} \, b^{2} {\left (\frac {5 \, e^{\left (-3 \, d x - 3 \, c\right )}}{d {\left (5 \, e^{\left (-2 \, d x - 2 \, c\right )} + 10 \, e^{\left (-4 \, d x - 4 \, c\right )} + 10 \, e^{\left (-6 \, d x - 6 \, c\right )} + 5 \, e^{\left (-8 \, d x - 8 \, c\right )} + e^{\left (-10 \, d x - 10 \, c\right )} + 1\right )}} - \frac {2 \, e^{\left (-5 \, d x - 5 \, c\right )}}{d {\left (5 \, e^{\left (-2 \, d x - 2 \, c\right )} + 10 \, e^{\left (-4 \, d x - 4 \, c\right )} + 10 \, e^{\left (-6 \, d x - 6 \, c\right )} + 5 \, e^{\left (-8 \, d x - 8 \, c\right )} + e^{\left (-10 \, d x - 10 \, c\right )} + 1\right )}} + \frac {5 \, e^{\left (-7 \, d x - 7 \, c\right )}}{d {\left (5 \, e^{\left (-2 \, d x - 2 \, c\right )} + 10 \, e^{\left (-4 \, d x - 4 \, c\right )} + 10 \, e^{\left (-6 \, d x - 6 \, c\right )} + 5 \, e^{\left (-8 \, d x - 8 \, c\right )} + e^{\left (-10 \, d x - 10 \, c\right )} + 1\right )}}\right )}
\]
[In]
integrate(csch(d*x+c)^3*(a+b*tanh(d*x+c)^3)^2,x, algorithm="maxima")
[Out]
-2*a*b*(arctan(e^(-d*x - c))/d - (e^(-d*x - c) - e^(-3*d*x - 3*c))/(d*(2*e^(-2*d*x - 2*c) + e^(-4*d*x - 4*c) +
1))) + 1/2*a^2*(log(e^(-d*x - c) + 1)/d - log(e^(-d*x - c) - 1)/d + 2*(e^(-d*x - c) + e^(-3*d*x - 3*c))/(d*(2
*e^(-2*d*x - 2*c) - e^(-4*d*x - 4*c) - 1))) - 8/15*b^2*(5*e^(-3*d*x - 3*c)/(d*(5*e^(-2*d*x - 2*c) + 10*e^(-4*d
*x - 4*c) + 10*e^(-6*d*x - 6*c) + 5*e^(-8*d*x - 8*c) + e^(-10*d*x - 10*c) + 1)) - 2*e^(-5*d*x - 5*c)/(d*(5*e^(
-2*d*x - 2*c) + 10*e^(-4*d*x - 4*c) + 10*e^(-6*d*x - 6*c) + 5*e^(-8*d*x - 8*c) + e^(-10*d*x - 10*c) + 1)) + 5*
e^(-7*d*x - 7*c)/(d*(5*e^(-2*d*x - 2*c) + 10*e^(-4*d*x - 4*c) + 10*e^(-6*d*x - 6*c) + 5*e^(-8*d*x - 8*c) + e^(
-10*d*x - 10*c) + 1)))
Giac [A] (verification not implemented)
none
Time = 0.39 (sec) , antiderivative size = 192, normalized size of antiderivative = 1.79
\[
\int \text {csch}^3(c+d x) \left (a+b \tanh ^3(c+d x)\right )^2 \, dx=\frac {60 \, a b \arctan \left (e^{\left (d x + c\right )}\right ) + 15 \, a^{2} \log \left (e^{\left (d x + c\right )} + 1\right ) - 15 \, a^{2} \log \left ({\left | e^{\left (d x + c\right )} - 1 \right |}\right ) - \frac {30 \, {\left (a^{2} e^{\left (3 \, d x + 3 \, c\right )} + a^{2} e^{\left (d x + c\right )}\right )}}{{\left (e^{\left (2 \, d x + 2 \, c\right )} - 1\right )}^{2}} + \frac {4 \, {\left (15 \, a b e^{\left (9 \, d x + 9 \, c\right )} + 30 \, a b e^{\left (7 \, d x + 7 \, c\right )} - 20 \, b^{2} e^{\left (7 \, d x + 7 \, c\right )} + 8 \, b^{2} e^{\left (5 \, d x + 5 \, c\right )} - 30 \, a b e^{\left (3 \, d x + 3 \, c\right )} - 20 \, b^{2} e^{\left (3 \, d x + 3 \, c\right )} - 15 \, a b e^{\left (d x + c\right )}\right )}}{{\left (e^{\left (2 \, d x + 2 \, c\right )} + 1\right )}^{5}}}{30 \, d}
\]
[In]
integrate(csch(d*x+c)^3*(a+b*tanh(d*x+c)^3)^2,x, algorithm="giac")
[Out]
1/30*(60*a*b*arctan(e^(d*x + c)) + 15*a^2*log(e^(d*x + c) + 1) - 15*a^2*log(abs(e^(d*x + c) - 1)) - 30*(a^2*e^
(3*d*x + 3*c) + a^2*e^(d*x + c))/(e^(2*d*x + 2*c) - 1)^2 + 4*(15*a*b*e^(9*d*x + 9*c) + 30*a*b*e^(7*d*x + 7*c)
- 20*b^2*e^(7*d*x + 7*c) + 8*b^2*e^(5*d*x + 5*c) - 30*a*b*e^(3*d*x + 3*c) - 20*b^2*e^(3*d*x + 3*c) - 15*a*b*e^
(d*x + c))/(e^(2*d*x + 2*c) + 1)^5)/d
Mupad [B] (verification not implemented)
Time = 3.77 (sec) , antiderivative size = 561, normalized size of antiderivative = 5.24
\[
\int \text {csch}^3(c+d x) \left (a+b \tanh ^3(c+d x)\right )^2 \, dx=\frac {a^2\,{\mathrm {e}}^{c+d\,x}}{d-d\,{\mathrm {e}}^{2\,c+2\,d\,x}}+\frac {136\,b^2\,{\mathrm {e}}^{c+d\,x}}{15\,\left (d+3\,d\,{\mathrm {e}}^{2\,c+2\,d\,x}+3\,d\,{\mathrm {e}}^{4\,c+4\,d\,x}+d\,{\mathrm {e}}^{6\,c+6\,d\,x}\right )}+\frac {32\,b^2\,{\mathrm {e}}^{c+d\,x}}{5\,\left (d+5\,d\,{\mathrm {e}}^{2\,c+2\,d\,x}+10\,d\,{\mathrm {e}}^{4\,c+4\,d\,x}+10\,d\,{\mathrm {e}}^{6\,c+6\,d\,x}+5\,d\,{\mathrm {e}}^{8\,c+8\,d\,x}+d\,{\mathrm {e}}^{10\,c+10\,d\,x}\right )}-\frac {a^2\,\ln \left (4\,a^6\,{\mathrm {e}}^{d\,x}\,{\mathrm {e}}^c-16\,a^4\,b^2-4\,a^6+16\,a^4\,b^2\,{\mathrm {e}}^{d\,x}\,{\mathrm {e}}^c\right )}{2\,d}+\frac {a^2\,\ln \left (4\,a^6+16\,a^4\,b^2+4\,a^6\,{\mathrm {e}}^{d\,x}\,{\mathrm {e}}^c+16\,a^4\,b^2\,{\mathrm {e}}^{d\,x}\,{\mathrm {e}}^c\right )}{2\,d}-\frac {2\,a^2\,{\mathrm {e}}^{c+d\,x}}{d-2\,d\,{\mathrm {e}}^{2\,c+2\,d\,x}+d\,{\mathrm {e}}^{4\,c+4\,d\,x}}-\frac {8\,b^2\,{\mathrm {e}}^{c+d\,x}}{3\,\left (d+2\,d\,{\mathrm {e}}^{2\,c+2\,d\,x}+d\,{\mathrm {e}}^{4\,c+4\,d\,x}\right )}-\frac {64\,b^2\,{\mathrm {e}}^{c+d\,x}}{5\,\left (d+4\,d\,{\mathrm {e}}^{2\,c+2\,d\,x}+6\,d\,{\mathrm {e}}^{4\,c+4\,d\,x}+4\,d\,{\mathrm {e}}^{6\,c+6\,d\,x}+d\,{\mathrm {e}}^{8\,c+8\,d\,x}\right )}+\frac {2\,a\,b\,{\mathrm {e}}^{c+d\,x}}{d+d\,{\mathrm {e}}^{2\,c+2\,d\,x}}-\frac {4\,a\,b\,{\mathrm {e}}^{c+d\,x}}{d+2\,d\,{\mathrm {e}}^{2\,c+2\,d\,x}+d\,{\mathrm {e}}^{4\,c+4\,d\,x}}-\frac {a\,b\,\left (\ln \left (32\,a^3\,b^3\,{\mathrm {e}}^{d\,x}\,{\mathrm {e}}^c+8\,a^5\,b\,{\mathrm {e}}^{d\,x}\,{\mathrm {e}}^c-a^5\,b\,8{}\mathrm {i}-a^3\,b^3\,32{}\mathrm {i}\right )\,1{}\mathrm {i}-\ln \left (32\,a^3\,b^3\,{\mathrm {e}}^{d\,x}\,{\mathrm {e}}^c+8\,a^5\,b\,{\mathrm {e}}^{d\,x}\,{\mathrm {e}}^c+a^5\,b\,8{}\mathrm {i}+a^3\,b^3\,32{}\mathrm {i}\right )\,1{}\mathrm {i}\right )}{d}
\]
[In]
int((a + b*tanh(c + d*x)^3)^2/sinh(c + d*x)^3,x)
[Out]
(a^2*exp(c + d*x))/(d - d*exp(2*c + 2*d*x)) + (136*b^2*exp(c + d*x))/(15*(d + 3*d*exp(2*c + 2*d*x) + 3*d*exp(4
*c + 4*d*x) + d*exp(6*c + 6*d*x))) + (32*b^2*exp(c + d*x))/(5*(d + 5*d*exp(2*c + 2*d*x) + 10*d*exp(4*c + 4*d*x
) + 10*d*exp(6*c + 6*d*x) + 5*d*exp(8*c + 8*d*x) + d*exp(10*c + 10*d*x))) - (a^2*log(4*a^6*exp(d*x)*exp(c) - 1
6*a^4*b^2 - 4*a^6 + 16*a^4*b^2*exp(d*x)*exp(c)))/(2*d) + (a^2*log(4*a^6 + 16*a^4*b^2 + 4*a^6*exp(d*x)*exp(c) +
16*a^4*b^2*exp(d*x)*exp(c)))/(2*d) - (2*a^2*exp(c + d*x))/(d - 2*d*exp(2*c + 2*d*x) + d*exp(4*c + 4*d*x)) - (
8*b^2*exp(c + d*x))/(3*(d + 2*d*exp(2*c + 2*d*x) + d*exp(4*c + 4*d*x))) - (64*b^2*exp(c + d*x))/(5*(d + 4*d*ex
p(2*c + 2*d*x) + 6*d*exp(4*c + 4*d*x) + 4*d*exp(6*c + 6*d*x) + d*exp(8*c + 8*d*x))) + (2*a*b*exp(c + d*x))/(d
+ d*exp(2*c + 2*d*x)) - (4*a*b*exp(c + d*x))/(d + 2*d*exp(2*c + 2*d*x) + d*exp(4*c + 4*d*x)) - (a*b*(log(32*a^
3*b^3*exp(d*x)*exp(c) - a^3*b^3*32i - a^5*b*8i + 8*a^5*b*exp(d*x)*exp(c))*1i - log(a^5*b*8i + a^3*b^3*32i + 32
*a^3*b^3*exp(d*x)*exp(c) + 8*a^5*b*exp(d*x)*exp(c))*1i))/d